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We revisit the deduction of the exit probability of the one dimensional Sznajd model through the Kirkwood 
approximation [F. Slanina et al, Europhys. Lett. 82, 18006 (2008)]. This approximation is peculiar in that in 
spite of the agreement with simulation results [F. Slanina et ah, Europhys. Lett. 82, 18006 (2008); R. Lambiotte 
and S. Redner, Europhys. Lett. 82, 18007 (2008); A. M. Timpanaro and C. P. C. Prado, 89, 052808 (2014)] 
the hypothesis about the correlation lenghts behind it are inconsistent and fixing these inconsistencies leads to 
the same results as a simple mean field. We use an extended version of the Sznajd model to test the Kirkwood 
approximation in a wider context. This model includes the voter, Sznajd and “United we Stand, Divided we 
Fall” (USDF) models [R. A. Holley and T. M. Liggett, Ann. Prob. 3, 643 (1975); K. Sznajd-Weron and J. 
Sznajd, Int. Joum. Mod. Phys. C 11, 1157 (2000)] as different parameter combinations, meaning that some 
analytical results from these models can be used to evaluate the performance of the Kirkwood approximation. 
We also compare the predicted exit probability with simulation results for networks with 10^ sites. The results 
show clearly the regions in parameter space where the approximation gives accurate predictions, as well as 
where it starts failing, leading to a better understanding of its reliability. 


1. INTRODUCTION 

In this work we investigate the reliability of the Kirkwood 
approximation, as used in fl] as a tool to study agent models 
defined in one dimensional lattices. This approximation is a 
type of mean field approach that keeps pair correlations but 
makes some unusual h^othesis about how the correlations 
decay. It was used in (^l to study the exit probability of the 
Sznajd opinion propagation model (the probability of reach¬ 
ing one of the two possible absorbing states as a function of 
the initial conditions) and a controversy about wether these 
results are indeed correct emerged in a series of papers |[ll- 
1^. Simulation results later stablished that the predicted exit 
probability is impressively acurate (it has been confirmed for 
simulations in networks with up to 3.16 x 10^ sites 1@]) and the 
same expression has since been obtained by a non mean field 
treatment that makes some hypothesis about how large groms 
of sites having the same opinion interact with each other 
Also, the fact that a mean field treatment manages to give an 
acurate answer in an one dimensional problem raised some 
speculations about wether there is something special about the 
way correlations behave in this model or if it is some kind of 
coincidence 101 ■ 

Motivated by this, we study the Kirkwood approximation 
in a wider context. We will focus on using this approximation 
to study the exit probability of a generalization of the Sznajd 
model that includes also the voter and “United we Stand, Di¬ 
vided we Fall” (USDF) models llsl-fl^. the reasoning being 
that this would create a parameter space where we can study 
the model. As we know that the approximation is acurate for 
some parameter values and expect this to be true for small 
perturbations of these parameters, this would lead to regions 
where the approximation works and regions where it fails, im¬ 
proving the understanding of what is behind it. 
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We will first define precisely the model used and make cal¬ 
culations for the exit probability (sectionlHli and later compare 
the obtained expression with simulation results (sectionlIII]i. 

A. The approximation in the Sznajd model 

The Sznajd model is an opinion propagation model, derived 
from the Ising model. In the situations that we will consider, 
the model is defined in an one dimensional chain with N sites 
and periodic boundary conditions. Each of its sites represents 
a person that can have one of two possible spin states, repre¬ 
senting opinions (denoted + and —). The time evolution is 
defined as follows; 

• At each time step, choose at random a pair of neigh¬ 
bouring sites i and i + 1. 

• If their opinions ai and ai+i are equal then we choose 
one of the neighbours of the pair, i — 1 or i + 2 and 
change its opinion to the pair opinion ai. 

• Otherwise, if ai ^ ai+i, nothing happens and we move 
on to the next time step. 

The Kirkwood approximation can be summed up as 

• Break correlations by pairing neighbouring spins and 
considering spins in a pair independent from spins in 
other pairs: 

{aia2a3) ~ {aia2){a3) 

• However, when we have only 2 spins we consider then 
correlated as if they were nearest neighbours: 

(cricr4) ~ {aia2). 

This unusual aspect of considering correlation lenghts to be 
short and long at the same time, depending on which corre¬ 
lation is being considered casts doubt about wheter the ap¬ 
proximation scheme even makes sense. However as we will 
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show, assuming a more coherent way of reducing the correla¬ 
tions (namely ((Ti(Ti+„) ~ {ai){ai+n)) reduces the results to 
the same as a simple mean field approximation (that gives a 
simpler, but inacurate, exit probability). 

We start by considering the master equation of the model. 
Since the only changes that are possible in each time step are 
single spin flips, then 


Cl = (CTi), C2{k) = (a^at+k), 

C 3 = (criCri+icri+ 2 ) and C 4 = ((Ti(Ti+icri+ 2 (Ti+ 3 ) ( 6 ) 
it follows from equations 0 ] and |5] that 


P{a) = Y, {W,{a^)P{a^) - W,{a)P{a)) , (1) 

i 

where a denotes the state of the chain, ct* is the state a flipping 
the site i, is the probability of going from state a to 

state cr® in one time step and the sum over i runs through all the 
sites in the chain. Consider then a function ip{a). Its ensemble 
average is given by 


(7 G 

Substituting the master equation and reordering the sums 
leads to 


= X X = 

G i 




( 2 ) 


We will be particularly interested in studying the average of 
spin products 


fc=i 

As al = ai(l — 2dij), substituting eq[3]in|2]leads to 


r Cl = Cl - C3 

1 C2(1) = 1-C2(1)+C2(3)-C4 

Applying the Kirkwood approximation leads us to the follow¬ 
ing 


C3~C2(l)Ci 

C^c^C2{l)^ (8) 

C2(3) ~C2(1) 


Since C 2 (l) is the only C 2 (k) that remains we will shorthand 
it to C 2 . It follows that 


r Cl = Ci(i-c2) 

I C2 = 1 - c| 


(9) 


We want to calculate the exit probability, which is the prob¬ 
ability of reaching the absorbing state where all spins equal 
-f as a function of the initial proportion p of sites -f (assum¬ 
ing also that the initial condition is completely uncorrelated). 
Since Ci = 1 in the all -f state and Ci = — 1 in the all — state 
(which are the only relevant absorbing states), we can obtain 
the exit probability by calculating the limit 


lim 

t—^OO 


1 -f Ci(f) 
2 


and using Ci(t = 0) = 2p — 1 and C 2 (t = 0) = Ci(f = 
0)^ = (2p — 1)^ as initial conditions. Doing this leads us to 
the following expression 


\ —1 / k —1 \ —1 / 

So to obtain the time evolution equations we must And Wi and 
use equation 0 ] to obtain an hierarchy of equations that will 
be truncated by the Kirkwood approximation. For the Sznajd 
model we have 


- crii.o-i+1 + CTi-i) - cri{a^+2 + cri-2)+ 

o 

+ 0'i+l<^i+2 + O'i-lO'i-l)- (5) 


Defining the following spin averages 


E{p) 


p2 + (l-p)2’ 


( 10 ) 


which matches the simulation results fairly well. However, as 
we pointed already, the Kirkwood approximation is inconsis¬ 
tent in the way it reduces correlations and this inconsistency is 
crucial for the result in eq[T 0 ] Obviously, if we assume that the 
correlations have longer ranges than just the first neighbour 
then we would need to treat C 2 (3), Cs and Ca as separate vari¬ 
ables and derive more equations through relation 01 involving 
more complicated spin averages and yielding a much more 
complex treatment. On the other hand if we reduce C 2 {k) in 
a way that is coherent with how the other correlations were 
treated then one should use C 2 (k) ~ Cf when /c > 1. If we 
substitute ( 72 ( 3 ) = Cf in eqs| 2 ]we get this system of equa¬ 
tions instead 
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r Cl = Ci(i-c2) 

\ C72 = 1 + C2 - C 2 - Cl 


( 11 ) 


and integrating these new equations leads to an exit probabil¬ 
ity given by a step function 


for these spin variables 6a,b,c = i'^+ab+ac+bc)/^ = A(a, b, c), 
meaning that 


^i{a,b, c) 


A{aai,b(Ji_i,cai-2) -f A{aai,bai+i,cai+2) 
2 


can be used to match the different spin patterns (±, ±, ±), 
(T, ±, ±), (±, T, ±) and (=F, T, ±), following that 


£;(p) = 0(^p-0, 

which is the same result as using a simple mean field (Ci = 
Ci-Cf), as well as applying the Galam Unified Frame to this 
problem d. 


^ p{a,b,+)^^{a,b,+). (12) 

a , fe —±1 

Expanding eq [12] leads to 

Wi((T) = 2a-f/3 (Tj((Tj+i -fcTi-i) + ja^{ai+2 + 0 -^- 2 )+ 


II. THE EXTENSION OF THE MODEL 

The extension of the Sznajd model that we will use was 
proposed by Kondrat in |@] and is defined by the following 
evolution rules 

• We choose a site i at random and take either the two 
neighbours to the left or the two to the right (with 50% 
probability each) to form a triplet. We will denote the 
nearest neighbour j and the second nearest k. 

• We flip site i {ai —J- —Ui) with probability p(CTi, aj,ak) 
and move on to the next time step. 

It follows that the model is completely specified by the p 
function. Since we have 2 possible spin states for each site 
there are 8 possible triplets {cTi,aj,ak). However, we are only 
interested in the situation where the + and — spins behave in 
a symmetric way, so we impose the constraint 


+ S{a^+ia^+2 + cri-icri- 2 ) (13) 

where a, /3 ,7 and 6 are the following parameter combinations 

P1+P2+P3+ Pi 
^- 

^ ^ P1-P2-P3+ Pi 


P1-P2+P3- Pi 



. Pl+P 2 -P 3 -Pi 

-8-■ 

Repeating what we did for the Sznajd model leads to 


— p{ ^ k) 1 

leaving us with 4 parameters to characterize p. These param¬ 
eters can be given an interpretation in terms of the type of 
opinion change they represent: 


Cl = -4((a + ^ + 7 )Ci + , 5 C 3 ) 

C 2 (l) = -4(/3+(2a + 7 )C 2 (l) + (/3 + ( 5 )C 2 ( 2 ) 
+7C2(3)+(5C4) 

(14) 

and after applying the approximation 


Pi — F( + ; +; +); 
P 2 =P(-,+,+), 
P 3 =Pi + ,-,+), 
Pi =Pi-,-,+), 


contrarian behaviour 
conformity behaviour 
conformity behaviour 
disagreement propagation. 


Moreover, 3 known models correspond to specific parameter 
choices: 


P2 = i,Pi = P3 = Pi = 0 leads to the Sznajd model 

F 2 = P 3 = Ij Pi = P4 = 0 leads to the voter model 

P 2 = P 4 = IjPi = P 3 = 0 leads to the USDF model. 


f Cl — — 4 Ci(a-f/ 3 -f 7-f 5C2) 

\ C 2 =-4(/3 + (2a +/3 + 27 + 5 )C 2 + ^C|) ^ ^ 

In order for the model to have an exit probability we need it to 
order, meaning that C 2 must go to 1 as time goes by. However, 
substituting C 2 = 1 in the equation for its derivative yields 
C 2 = — 8 (a -I- /3 -f 7 -f b) = —4pi. Hence, we must restrict 
ourselves to models such that pi = 0 (which makes sense, 
because the ferromagnetic states are not absorbing otherwise). 
As this means that a + j3 + ^ + 5 = Q the equations can be 
further simplified: 


We want to And a system of equations similar to eqs|2 mean¬ 
ing we must obtain an expression for Wi (a). We first note that 


r Cl =45 Ci(1-C2) 

\ C 2 = -4(/3 - (/3 + S)C2 + SCI) 
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We also want that the solution (72 = 1 be attractive, which 
means 25 — {fi -\- 5) >0<^5>/3<^p2 > P 4 - With 
these two conditions in place the integration of the equations 
[T^ becomes very similar to what was done in the case of the 
Sznajd model and we can obtain the following exit probability 


E{p) 


P(1 - P)P3+ P'^{P2 - Pi) 

2p{l - p)p3 + (p2 + (1 - p)2)(p2 - Pi) ■ 


(17) 


Once again we can try substituting (72(2) and (72(3) by Cf 
and the exit probability obtained is the same as in the case 
of a mean field with no correlations (i.e. integrating Ci = 
4i5(7i(l — Cl)), with a step function for b < 0, a linear exit 
probability for b = 0 and the model failing to order if (5 > 0. 


III. SIMULATION RESULTS 

We now make simulations of the model to compare the re¬ 
sults with the expressionfm Since we are restricting ourselves 
to Pi =0 and p 2 > Pi we have 3 free parameters: p 2 , ps and 
Pi. However the parameter space we need to explore has only 
2 dimensions since the model with parameters p' = ppi dif¬ 
fers from the model with parameters pi only on the time scale 
where things happen, so that the limit for large times remains 
the same (this can be seen in equation[T2]that is invariant by a 
rescaling of the parameters). As a consequence, we only need 
to analyse two different cases: P 2 = 1 with pa and pi free, 
and p3 = 1 with p2 > Pi- 

Before making the simulations we take a closer look at 
the prediction given by the Kirkwood approximation, eq [17] 
Since the expression is invariant by a reescaling of the pa¬ 
rameters and only 2 parameter combinations appear (pa and 
P 2 — Pi) we can rewrite it using only one parameter tjj and the 
dependance of the exit probability with p would be the same 
along lines with constant tp. We choose 


P2 -Pi 

that varies from 0 when P 2 = 1 and pa = 0 and diverges in the 
limit p 2 Pi, when the Kirkwood approximation becomes 
problematic. The exit probability in eafTTlbecomes 


E{p) 


tppjl - p) + p'^ 

2^p(l - p) -I- p2 -I- (1 - p)2 


(18) 


and the curves ip = constant can be seen in figure [T] Note 
the singularity that arises for P 2 = P 4 = 1 and pa = 0, corres¬ 
ponding to the USDF model (that has antiferromagnetic states 
with a non-zero probability of being reached in addition to the 
ferromagnetic ones, meaning that the exit probability as we 
defined is not adequate in this case). 

We also know that the model with parameters such that 
P 2 = Pd. E Pi (corresponding to ip = 1) must have a linear exit 
probability (excluding the USDF case). This happens because 
the model can be rewritten (up to a time scale) as follows: 



FIG. 1: Curves along which ip = is a constant in 

parameter space. Along each curve the functional form of 
E{p) predicted by the Kirkwood approximation is the same. 
Note the singularity for P 2 = P 4 = 1 and pa = 0 that corres¬ 
ponds to the USDF model. 


• At each time step choose a site i. 

• With probability P^/p^, i copies the opinion of one of its 
first neighbours. 

• Otherwise (hence with probability Pi/p 2 ) i copies the 
opinion of one of its second neighbours. 

Since both of the possible interactions conserve the magneti¬ 
zation (N+ — N-) on average, the ensemble average of p after 
very long times (corresponding to the exit probability in these 
cases) is the same as in the initial condition, implying E = p. 
Because of this we didn’t made simulations in this case, using 
the analytical result instead. 

For the simulations we used linear chains with periodic 
conditions and 10^ sites. While this is a relatively small 
network size it is enough to get a rough picture of how the 
exit probability behaves (the main reason larger network sizes 
were not used is that the techniques used to speed up the 
results for the Sznajd model (@] are not as efficient when 
Pi ^ 0). For each parameter set we made simulations using 
p = 0.1, 0.2, 0.3, 0.4 and 0.45, calculating the absolute dif¬ 
ference between the simulation results for the exit probability 
and the one predicted with the Kirkwood approximation and 
to visualize these in two dimensions we made the average of 
this absolute difference (that works as an estimate of the ab¬ 
solute error). Finaly, depending on the parameters, we made 
between 10 ® and 10 ® simulations (depending mainly on how 
long the simulations took and how much the measured exit 
probability deviated from the one predicted by the Kirkwood 
approximation). The results obtained are summed up in figure 

El 

We can distinguish easily that the Kirkwood approximation 
starts failing as pi P 2 , as well as two distinct regions where 
the approximation has greater accuracy. One of them is the 
small region around the Sznajd model {p 2 = l,pa = Pi = 0, 
corresponding to the lower left corner in the graph of figure|2|i 
and the other contains the vicinity of the voter model {p 2 = 
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FIG. 2; The graph shows the estimated absolute error of the 
Kirkwood approximation for each of the simulated parameter 
sets. Brighter colors indicate smaller differences, while darker 
colors indicate larger differences (values larger than 0.007 are 
all in black). The marked points indicate the data we have 
and that was interpolated to yield the graph, while the hashed 
region corresponds to regions where we had no data. The line 
P 2 = P 3 + Pi is also marked, where the error was set to 0. 
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FIG. 3: The graph shows the estimated absolute error of equa¬ 
tion [18] with a fitted value of 'ip, for each of the simulated 
parameter sets. Brighter colors indicate smaller differences, 
while darker colors indicate larger differences (values larger 
than 0.007 are all in black). The marked points indicate the 
data we have and that was interpolated to yield the graph, 
while the hashed region corresponds to regions where we had 
no data. The line p 2 = P 3 +P 4 is also marked, where the error 
was set to 0. 



FIG. 4: Simulation results (*) and the Kirkwood approxima¬ 
tion (darker curve) for the parameter choice p 2 = Ps = 1, 
P 4 = 0.875. The curve E = pis drawn for reference (lighter 
line). Statistical errors are too small to see. 


half of the parameter space), where the minority proportion 
always increases, but the model still orders. It is important 
to note that this means that an uncorrelated initial condition 
will see an increase of the minority opinion, while this ceases 
to be true once correlations form (allowing the system to or¬ 
der). This must be an effect of simulating the model in one 
dimension and shouldn’t happen in networks offering more 
possibilities for the opinions to mix. Figures|4|and|5]show the 
Kirkwood approximation and simulation results for two pa¬ 
rameter choices (one where the approximation is acurate and 
other where it is not). 


IV. CONCLUSION 


P 3 = l,pi = 0) as well as the line P 2 = P 3 + Pi and part of 
P3 =P2 +Pi- 

Nevertheless, equation [18] alone gives a good description 
of the exit probability in most of the parameter space if we 
allow 'ip to take values different than the ones predicted by the 
Kirkwood approximation. We repeated the comparison made 
in figure [ 2 ] but now using a fitted value of 'ip instead of the 
one predicted by the approximation. This comparison can be 
found in figure[3] where we can see that most of the parameter 
space seems to obey the functional form in eq[T8]even though 
only a part of it matches the exact prediction of the Kirkwood 
approximation. 

Finally we call to attention the peculiar aspect of the exit 
probability when 'ip > 1 (which corresponds to more than 


In this work we approached the controversy about the exit 
probability of the Sznajd model through a different angle, 
aiming at the validity of the Kirkwood approximation itself. 
Our simulations show some clear regimes where the simula¬ 
tion does not yield accurate results, but nevertheless it seems 
to give correct qualitative results (like the region in the pa¬ 
rameter space where the system orders and the general func¬ 
tional form of the exit probability), as would be expected from 
a mean field treatment. We have also shown that the way cor¬ 
relations are approximated appears to be inconsistent and that 
if we make instead consistent choices, the conclusions about 
the exit probability and the parameters for which the system 
orders match the ones obtained by a simple mean field approx¬ 
imation. On the other hand the approximation also gives re- 
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markably acurate predictions in some parameter regions. This 
acuracy was expected for the line P 2 = P 3 + Pi as we pointed 
out, because the magnetization is conserved in average so the 
correlations are not all that important. The large region in the 
center of figure |2] however was completely unexpected and 
may give a hint about when the inconsistencies in the treat¬ 
ment of the correlations turn out to be really relevant. In par¬ 
ticular, part of the region where p 2 < Pa -b Pi has a peculiar 
shape for the exit probability curve that has been confirmed by 
the simulations (see hgure [S]) and that would imply that mi¬ 
nority opinions become represented disproportionately more 
than majority opinions. Situations like this should not intu¬ 
itively happen unless the model failed to order. In our opin¬ 
ion this is an effect of simulating the model in one dimension, 
which severely curtails how much the opinions can mix. What 
is also remarkable is that the Sznajd model seems to be in a 
small region of acuracy,which could indicate that there really 
is something special about the way correlations behave in it. 


FIG. 5; Simulation results (*) and the Kirkwood approxi¬ 
mation (darker curve) for the parameter choice p 2 = 0.75, 
P 3 = 1, p 4 = 0.25. The curve E = p is drawn for reference 
(lighter line). Statistical errors are too small to see. 
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